
function Hv=H(alpha,beta)

k=[cos(alpha)*cos(beta);
   cos(alpha)*sin(beta);
   sin(alpha)];

r{1}=[0, 0, 0];
r{2}=[1/4, 0, 0];
r{3}=[1/8, -1/4, 0];
r{4}=[-1/8, 1/4, 0];
r{5}=[-1/4, 0, 0];
r{6}=[-1/8, 1/4, 0];
r{7}=[1/8, 1/4, 0];

phi=nan(7,1);

for i=1:7
    phi(i)=2*pi*r{i}*k;
end

j=sqrt(-1);
Hv=exp(j*phi);
end